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C \ Abstract 



The theory of interactions between lasers and cold trapped ions as it pertains to the 
design of Cirac-Zoller quantum computers is discussed. The mean positions of the 
trapped ions, the eigenvalues and eigenmodes of the ions' oscillations, the magnitude 
of the Rabi frequencies for both allowed and forbidden internal transitions of the ions 
and the validity criterion for the required Hamiltonian are calculated. Energy level 
data for a variety of ion species is also presented. 

PACS numbers: 32.80.Qk, 42.50.Vk, 89.80.+h 



1 Introduction 

A quantum computer is a device in which data can be stored in a network of quan- 
tum mechanical two-level systems, such as spin-1/2 particles or two level atoms. The 
quantum mechanical nature of such systems allows the possibility of a powerful new 
feature to be incorporated into data processing, namely the capability of performing 
logical operations upon quantum mechanical superpositions of numbers. Thus in a 
conventional digital computer each data register is, throughout any computation, al- 
ways in a definite state "1" or "0" ; however in a quantum computer, if such a device 
can be realized, each data register (or "qubit") will be in an undetermined quantum 
superposition of two states |1) and |0). Calculations would then be performed by ex- 
ternal interactions with the various two-level systems that constitute the device, in 
such a way that conditional gate operations involving two or more different qubits 
can be realized. The final result would be obtained by measurement of the quantum 
mechanical probability amplitudes at the conclusion of the calculation. Much of the 
recent interest in practical quantum computing has been stimulated by the discovery 
of a quantum algorithm which allows the determination of the prime factors of large 
composite numbers efficiently [|1|, and of coding schemes that, provided operations on 
the qubits can be performed within a certain threshold degree of accuracy, will allow 
arbitarily complicated quantum computations to be performed reliably regardless of 
operational error p|. 

So far, the most promising hardware proposed for implementation of such a de- 
vice seems to be the cold-trapped ion system devised by Cirac and ZoUer 0. Their 
design, which is shown schematically in figure 1, consists of a string of ions stored in a 
linear radio-frequency trap and cooled sufficiently that their motion, which is coupled 
together due to the Coulomb force between them, is quantum mechanical in nature. 
Each qubit would be formed by two internal levels of each ion, a laser being used 
to perform manipulations of the quantum mechanical probability amplitudes of the 
states; conditional two-qubit logic gates being realized with aid of the excitation or 
de-excitation of quanta of the ions' collective motion. For a more detailed description 
of the concept of cold-trapped ion quantum computation, the reader is referred to the 
article by Steane ^. 

There are two distinct possibilities for the choice of the internal levels of the ion: 
firstly the two states could be the ground state, and a metastable excited state of 
the ion (or more precisely, sublevels of these states); secondly the two states could be 
two nearly degenerate sublevels of the ground state. In the first case a single laser 



would suffice to perform the required operations; in the second, two lasers would be 
required to perform Raman transitions between the states, via a third level. Both of 
these schemes have advantages: the ffist, which I will refer to as the "single photon" 
scheme, has the great advantage of conceptual and experimental simphcity; the second, 
the "Raman scheme" , offers the advantages of a very low rate for spontaneous decay 
between the two nearly degenerate states and resilience against fluctuations of the 
phase of the laser. This later scheme was recently used by the group headed by Dr. D. 
J. Wineland at the National Institute of Science and Technology at Boulder, Colorado 
to realize a quantum logic gate using a single trapped Beryllium ion 0. 

In this paper I will discuss the theory of laser interactions with cold trapped ions 
as it pertains to the design of a Cirac-ZoUer quantum computer. I will conscentrate 
on the "single photon scheme" as originally proposed by those authors, although much 
of the analysis is also relevant to the "Raman scheme" . Fuller accounts of aspects of 
this are availible in the literature: see for example [^, 0, 0]; however the derivation 
of several results are presented here for the first time. I will also present relevant data 
gleaned from various sources on some species of ion suitable for use in a quantum 
computation. 

2 Equilibrium positions of ions in a linear trap 

Let us consider a chain of A^ ions in a trap. The ions are assumed to be strongly 
bound in the y and z directions but weakly bound in an harmonic potential in the x 
direction. The position of the m — th ion, where the ions are numbered from left to 
right, will be denoted Xm{t)- The motion of each ion will be influenced by an overall 
harmonic potential due to the trap electrodes, and by the Coulomb force exerted by all 
of the other ions. Hence the potential energy of the ion chain is given by the following 
expression: 

V = E ^Mu^^Utf + E ^ 1 ,,, ,.| , (2.1) 

^1 2 „^tli 87reo |a;„(t) - Xm{t)\ 

where M is the mass of each ion, e is the electron charge, Z is the degree of ionization 
of the ions, eo is the permitivity of free space and z/ is the trap frequency, which 
characterizes the strength of the trapping potential in the axial direction. Note that 
this is an unconventional use of the symbol z/, which often denotes frequency rather 
than angular frequency; following Cirac and Zoller, I will use u to denote the angular 
frequencies of the laser or the transitions between internal states of the ions, and z/ to 



denote angular frequencies associated with the motion of the ions. 

Assume that the ions are sufficiently cold that the position of the m — th ion can 
be approximated by the formula 

xUt) ~ x^^^ + qm{t) (2.2) 

where x^^ is the equilibrium position of the ion, and qm{t) is a small displacement. 
The equilibrium positions will be determined by the following equation: 

dV 



dXr 



= (2.3) 

(0) 



If we define the length scale i by the formula 



47reoMz/2 ' ^^'^^ 

and the dimensionless equilibrium position Um = x^m I ^i then (|2.3| ) may be rewritten 
as the following set of A^ coupled algebraic equations for the values of Um'- 

m-l 1 N 1 

Um-Y.1 ^^+ E 7 ^^ = {rn = l,2,..N). (2.5) 

For N = 2 and A^ = 3 these equations may be solved analytically: 

A^ = 2: m = -{l/2f\ «2 = (l/2)2/^ (2.6) 

N = 3: ui = -(5/4)l/^ U2 = 0, m = (5/4)^/3 _ (2.7) 

For larger values of A^ it is necessary to solve for the values of Um numerically. The 
numerical values of the solutions to these equations for 2 to 10 ions is given in table 
1. Determining the solutions for larger numbers of ions is a straightforward but time 
consuming task. 

By inspection, the minimum value of the spacing between two adjacent ions occurs 
at the center of the ion chain. Compiling the numerical data for the minimum value of 
the separation for different numbers of trapped ions, we find that it obeys the following 
relation: 

-rmniN) ^ ^ . (2.8) 

This relation is illustrated in figure 2. Thus the minimum inter-ion spacing for different 
numbers of ions is given by the following formula: 

/ ZV y/' 2.018 

This relationship is important in determining the capabilities of cold trapped ion quan- 
tum computers M. 



3 Quantum Fluctuations of the ions 

This section discusses the equations of motion which describe the displacements of the 
ions from their equihbrium positions. Because of the Coulomb interactions between 
the ions, the displacements of different ions will be coupled together. The Lagrangian 
describing the motion is then 



M 



N 



N 



L = ^J2 (^'»)^ ~oJ2 (InQn 



m=l 



n,m=l 






(3.i: 



where the subscript denotes that the double partial derivative is evaluated at g„ = 
gm = 0. The partial derivatives may be calculated explicitly to give the following 
expression: 



2 



N 
m=l 



N 
n,m=l 



(3.2) 



where 



/ir 



N 



1 + 2E 
p=l 

-2 



p=l l^m "^pl 



if n = m. 



if n ^ m. 



(3.3) 



l^m ^n I 

Since the matrix Anm is real, symmetric and non-negative definite, its eigenvalues 
must be non-negative. The eigenvectors 6^^ {p = 1,2, ■ ■ ■ A^) are therefore defined by 
the following formula: 



N 



EA„^6(f)=^,6(^)(p=l,...,iV) 



(3.4) 



n=l 



where fip > 0. The eigenvectors are assumed to be numbered in order of increasing 
eigenvalue and to be properly normalized so that 



N 

p=i 

N 
n=l 



(5r, 



Opq ■ 



(3.5) 
(3.6) 



The first eigenvector (i.e. the eigenvector with the smallest eigenvalue) can be 
shown to be 



b(i) 



{1,1,...,!}, /ii = 1 . 



(3.7) 



The next eigenvector can be shown to be 



b^^) = "72 {mi,W2,---,M7v}, /i2 = 3 . (3i 



N 
\m=l ) 



E^^ 



Higher eigenvectors must, in general, be determined numerically. Equations ( |3.6| ) and 
( |3.7| ) imply that 

N 

5:6(^) = ifp^l. (3.9) 

m=l 

For N = 2 and A^ = 3, the eigenvectors and eigenvalues may be determined 
algebraically: 

N = 2: b« = ^(1,1), /ii = l, 

(3.10) 
iV = 3: b« = ^(1>1,1)> /^i = l, 

b(2) = -^(-1,0,1), /i2 = 3, 

b(3) = i=(i,_2,l), /is = 29/5, (3.11) 

For larger values of A^, the eigenvalues and eigenvectors must be determined numeri- 
cally; their numerical values for 2 to 10 ions are given in table 2. 
The normal modes of the ion motion are defined by the formula 

Qpit) = E btkmit) . (3.12) 

m=l 

The first mode Qi{t) corresponds to all of the ions oscillating back and forth as if they 
were rigidly clamped together; this is referred to as the center of mass mode. The 
second mode Q2(t) corresponds to each ion oscillating with an amplitude proportional 
to its equilibrium distance form the trap center; This is called the breathing mode. The 
Lagrangian for the ion oscillations (|3.2| ) may be rewritten in terms of these normal 
modes as follows: 



2p^i 



where the angular frequency of the p — th mode is defined by 

^P = ^/Jh^■ (3-14) 

This expression imphes that the modes Qp are uncoupled. Thus the canonical mo- 
mentum conjugate to Qp is Pp = MQp and one can immediately write the Hamiltonian 

as 

H=^T.^P; + fT./rQl- (3-15) 

The quantum motion of the ions can now be considered by introducing the operators Q 



(op-aj) , (3.16) 



\ 2Mu. 
Pp^Pp = ^^^(ap + a\) . (3.17) 

where Qp and Pp obey the canonical commutation relation Qp, Pq = iMpq and the 
creation and annihilation operators d^ and dp obey the usual commutation relation 



dp, Clg 



Opq- 



Using this notation, the interaction picture operator for the displacement of the 
m — th ion from its equilibrium position is given by the formula: 

N 



p=l 



V 2MuN p^, 
where the coupling constant is defined by 



J2 s^ZHapC-^"''' - ale'^'n , (3-18 



_ _rNbiP) 
For the center of mass mode 



st^ = ^^^T^ ■ (3.19) 



1/4 
^P 



.(1) 



s\^^ = l v^ = v, (3.20) 

and for the breathing mode 



^m = ^— ^ TV^"™ z/2 = y3z/. (3.21) 



E 



u^ 



m 
\'m=l / 



^ There is some arbitarinees in the definition of the operators Pp and Qp, which is related to the 
arbitrariness of the phase of the Fock states. I have used the definitions given by Kittel (ref . |9[| , p. 16), 
which differs from that given in other texts on quantum mechanics (see for example, ref. |1C| , p. 183 
or ref. HI p. 36). 



4 Laser-Ion interactions 

I will now consider the interaction of a laser field with the trapped ions. The the- 
ory must take into acount both the internal and vibrational degrees of freedom of the 
ions. I will consider two types of transition between internal ionic levels: the familiar 
electric-dipole allowed (El) transitions and dipole forbidden electric quadrupole (E2) 
transitions. Electric quadrupole transitions have been considered in detail by Freedhoff 
jl^,[jl3|. The reason for considering forbidden transitions is that they have very long 
decay lifetimes; spontaneous emission will destroy the coherence of a quantum com- 
puter, and therefore is a major limitation on the capabilities of such devices [Q,[^. 



Magnetic dipole (Ml) transitions, which also have long lifetimes, tend only to occur 
between sub-levels of a configuration, and will therefore require long wavelength lasers 
in order to excite them. As it is necessary to resolve individual ions in the trap using 
the laser, the use long wavelengths will seriously degrade performance. More highly 
forbidden transitions are also a possibility for use in a quantum computer. In partic- 
ular, there is a octupole allowed (E3) transition at 467nm with a theoretical lifetime 
of 1.325 X 10^ secfl^, which has recently been observed at the National Physical 
Laboratory at Teddington, England |]I6| . 



The interaction picture Hamitonians for electric dipole (El) and electric quad- 
rupole (-E2) transitions of the m-th ion, located at x^ are 

i^p) = teY.UMN\N){M\{N\f,\M)A,i^,^,t)e''^^'^\ (4.1) 

MN 

^(E2) ^ |^^^^|iV)(M|(iV|f,f,|M)a,A-(x™,t)e''--^-*, (4.2) 

^ MN 

where Aj{'x,t) is the j — th component of the vector potential of the laser field, di 
denotes differentiation along the i — th direction and summation over repeated indices 
(z, j = x,|/, z) is implied; fj is the i — th component of the position operator for the 
valence electron of the ion; {|A^)} is the set of all eigenstates of the unperturbed ion 
and the transition frequency is ujmn = ^m — ^n where the energy of the N — th state 
is huj]s[ . 

For a laser beam in a standing wave configuration (see fig. 1), propagating along 
a direction specified by the unit vector n, the vector potential and its derivative are 
given by the formulas 

E r - 1 
Ai(x^,t) = -ei— sm fcCm(t) e*^* + c.c, (4.3) 



OiAjyx.m, tj 



E 



—Uiej — cos 
c 



Km{t) 



e*"* + c.c. 



(4.4) 



In (|4.4| ), I have approximated the laser beam as a plane wave, e being the polarization 
vector, E is the amplitude of the electric field, uj is the laser frequency and k = uj/c is 
the wavenumber. The operator (m{t) is the distance between the m — th ion and the 
plane mirror used to form the standing wave. 

If we restrict our consideration to just two states, 1 1) and |2) , and make the rotating 
wave equation, the interaction Hamiltonians may be rewritten as follows: 



rim 



hO,r, sin 






f^smyt) 



J(tA~ 



'|1)(2| 



fcC„(t)le^(*^~^)|l)(2| 



+ 
+ 



h.a. , 
h.a. 



where the detuning is A 



UJ 



UJ21 and the Rabi frequencies are given by 



a 



(El) 



Q, 



(E2) 



eE 

T 

eEuj2i 



2hc 



{l\firj\2)einj 



(4.5) 
(4.6) 

(4.7) 
(4.8) 



If the standing wave of the laser is so contrived that the equilibrium position of 
the m — th ion is located at a node, i.e., the electric field strength is zero, then 



(rn{t) = IX + COsOqmit) 



(4.9) 



where / is some integer, A is the wavelength, 6 is the angle between the laser beam and 
the trap axis and we have assumed that the fiuctuations of the ions transverse to the 
trap axis are negligible. In this case the two Hamiltonians become: 



H'l 



(El) 



hni^^hcoseqm{t)e'^'^- 



>+i7r) 



1)(2| + h.a. 



H 



(B2) 



;if]f')e*[*'^-'^-('+i/2)-]|l)(2| + h.a. 



(4.10) 
(4.11) 



where we have neglected terms involving qm{tY- It is convenient to write the dis- 
placement of the ion in terms of the creation and annihilation operators of the phonon 
modes, viz.: 



kcosOq^it) 



V 



N 






-iUrjt 



P=l 



aJe-^S, 



(4.12) 



where rj = Jhk'^cos'^6 / 2Mz/ is called the Lamb-Dicke parameter. 

Similarly if the standing wave is arranged so that the ion is at an antinode, i.e. 

(2/ - 1)A 



cut) 



+ coseq^{t) 



(4.13) 



then the Hamiltonians are: 

Hf^^ = ;ifi[,^'^e*(*^-'^+'")|l)(2| + h.a. , (4.14) 

^P) = ;ifif')A;cos^g^(t)e*[*^-^-('+V2)-]|i)(2| + h.a. (4.15) 

Thus we have two basic types of Hamiltonian: 

Hv = hQoe'^''^"''^^\l){2\ + h.a. , (4.16) 

Hu = hVtokcoseq^{t)e'^'^-'t"^'^\l){2\ + h.a., (4.17) 

where VLq stands for either VLq or ^q . 

By changing the node to an antinode, by moving the reflecting mirror for example, 
we can switch from one type of Hamiltonian to the other. In the first case, the laser 
beam will only interact with internal degrees of freedom of the ion, while in the second 
case the collective motion of the ions will be affected as well. 

5 Evaluation of the Rabi frequencies 

We can relate the matrix elements appearing in the definitions of the Rabi frequencies 
to the Einstein A coefficients for the transitions. In order to do this we will rewrite the 
matrix elements in terms of the Racah tensors: 

(l|r.|2)6. = Y. (l|rC«|2)cS'')e., (5.1) 

(l|f,f,|2)6,n,. = Y^ (l|r2cf)|2)cSfe,n,, (5.2) 

q=-2 



where we have used the fact that e-n = 0. The vectors c*^'') and the second rank tensors 
q| may be calculated quite easily; explicit expressions are given in the appendix. If 
we assume LS coupling, the states |1) and |2) are specified by the angular momentum 
quantum numbers; thus we will use the notation |1) = \jmj) and |2) = \j'm'j) , where 
j is the total angular momentum quantum number and rrij the magnetic quantum 
number of the lower state and j' is the total angular momentum quantum number 
and m'- the magnetic quantum number of the upper state. Using the Wigner-Eckart 
theorem (ref. |T^, section 11.4), the matrix elements may be rewritten as 

(1I..I2).. = iHlrC'^im ±[ J^^ I ,^)c«"... (5.3) 

{l|nr,|2><,n, = (lllr^CPIp) T. ( J,^ I £, ) c'ff.n,. (5.4) 



-m q m 
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the terms containing six numbers in brackets being Wigner 3 — j symbols (ref. |T1 
section 5.1), and (l||r''C^''^||2) being the reduced matrix element. The Einstein A 
coefficients for the two levels are given by the expressions: 



^12 



A 



(E2) 



12 



^ E |(iK^|2) 

c«^?2 ^ |(l|r2cf |2) 



15 



q=-2 



(5.5) 
(5.6) 



Using the Wigner-Eckart theorem again, these expressions reduce to the following: 



^12 



a{E2) 
^12 



Acak 



12 



cak\2 



(i|kc«l|2)r E 



q=-l 



J 1 f 

-rrij q m'. 



15 



[l\\r^C^'^\\2)f J: ( ^ ^ ^', 



(5.7) 
(5.8) 



These coefficients are the rates for spontaneous decay from the upper level |1) to 
the lower level |2). A simpler expression for the total rate of spontaneous decay of |2) 
to all of the sublevels of the lower state may be found by summing these rates over all 



values of m^: 



^(Sl) _ Y- a(E1) 



H2 



AiE2) 
^12 



- E A2 

m=-j 
3 



E A^ 



m=-j 



4cakl2 
3(2j'+l) 

cakl2 
15(2/ + 1) 



;i||rC«| 



(l|k'C(2)||2) 



(5.9) 
(5.10) 



These decay rates, which are the same for all of the sublevels of the upper level, are 
the quantities usually quoted in data tables. Using ( ^I7| - [4.8|) , (|5.3[ - |5.4|) and (|5l9| - |5.10|) , 
we then obtain the following formula for the Rabi frequencies: 



fio 



eE A 



12 



hy/ca\ kf. 



-0", 



(5.11) 



12 



where 



a 



(El) 



'3(2j' 






a 



(E2) 



'15(2/ + 1) 



2 

E 

q=~2 



J 2 f 

-rrij q m'. 



^e, 



(g) 



(5.12) 
(5.13) 
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The values of these quantities will be dependent on the choice of states of ions 
used for the upper and lower levels, and upon the polarization and direction of the laser 
beam. As a specific example, we will assume that the ions are in a weak magnetic field, 
which serves to define the z-direction of quantization. Furthermore, we will assume that 
the lower level |1) is the rrij = —1/2 sublevel of a ^5*1/2 ground state, the nucleus having 
spin zero. The upper level for the dipole transition is a sublevel of a ^Pi/2 state, while 
for the quadrupole transition it is a sublevel of a ^-D3/2 state. 



(El) 



n. 



"O 



e\E\ 



^12 



(El) 



r\ 4^^' ^^-^^^ 



T- (5-15) 



h \ 2cak 

6 Validity of Cirac and Zoller's Hamiltonian 



Equations (|4.12| ) and ( [4.17|) give the following expression for the Hamiltonian for the 



case when the laser standing wave is so configured that it can excite the vibration 
modes of the ions: 



Hu = '-^ E ^^^ i^P^^'""' - ^e'^"') e^^*^"^"^|l)(2| + h.a., (6.1) 



p=i 

In their paper P], Cirac and ZoUer assumed that the laser can interact with only 
the center of mass mode of the ions' fluctuations. This interaction forms a vitally 
important element of their proposed method for implementing a quantum controlled 
not logic gate. Thus they used a Hamiltonian of the following form [c.f. ref.[0, eq. (1)] 

^(cz) ^ ^rmo ^^^^.i,,t _ ^te-it) g^{tA-<^.)|^^^2| + h.a. (6.2) 

This is an approximate form of ( |6.1D , in which all of the other "extraneous" phonon 
modes have been neglected. We will now investigate under what circumstances these 
modes may be ignored. 

We will assume that the wavefunction for a single ion interacting with the laser 
beam may be written as follows: 

N N 

\^it)) = aoit)\l)\vac) +boit)\2)\vac) + Y.aM^)\lp) + J2bpit)\2)\l,), (6.3) 

p=i p=i 
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where |1) and |2) are the energy eigenstates of the m — th ion's internal degrees of 
freedom, |lp) is the state of the ions' collective vibration in which the p — th mode 
has been excited by one quantum, and \vac) is the vibrational ground state. To avoid 
ambiguity, the ket for the ion's internal state appears first, the ket for the vibrational 
state second. 

The equation of motion for this wavefunction is 

th^^\^it))=Hu\^{t)). (6.4) 

Using (|6.1| ), and assuming that one cannot excite states with two phonons, one obtains 
the following equations: 



^oV 



N 



«o = ^E^^^/^pW (6-5) 

$0 - 



N 



-^ E 4^^«.(^) (6.6) 



p=l 

dp = -z(z/p - i^i)ap ^s^)/3o(t) (6.7) 

y JM 

/3p = -^{u, + u,)P,-^s^Z^aoit) (6.8) 



We have assumed that A = — z/i, so that the laser is tuned to the specific sideband 
resonance required to perform Cirac and Zoller's universal gate operation (ref. 0, eq. 
3), namely the two level transition |l)^|li) ^> \2)^\vac). 

Since |ao(^)|) l/3o(^)| < 1; "we can consider the following upper limits on the ampli- 
tudes of the states which include excitation of "extraneous" phonon modes (i.e. phonon 
modes other than the center of mass mode): 

\ap{t)\<\Ap{t)l \f3p{t)\ < \Bp{t)\ (6.9) 

where 

Ao + ^{u,-uMp = -^st^ (6.10) 



5o + z(z/p + z/i)5p = -^/l\ (6.11) 



A^ 



Solving these equations one finds that 



\AM < ^#^I4^^I (6.12) 

\BAt)\ < J^'l A st\ (6.13) 
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Thus the total probabihty that "extraneous" modes are excited has the following upper 
limit: 



Pe.t = j:K{t)\'+mt)\'<2(^^^' 



p=2 



N 



1 






.(P)^2 



p=2 



(6.14) 



where we have used the definition of the mode frequencies (|3.14|) and the fact that the 
eigenvalue for the center of mass mode is /ii = 1. This quantity will be different for 
each ion in the string; taking its average value, we find 



N 



ext 



TiEP. 



N 



ext 



m=l 



< ^ ffi) ^"^'' 



(6.15) 



where we have used the definition of the coupling constants, ( |3.19|) and the orthonor- 
mality of the eigenvectors ( |3.6| ). The function S(A^) is defined by the formula 



N 



s(iv) = E 



/ip + i 



p=2 



(/ip- ly^' 



(6.16) 



This must be evaluated numerically by solving for the eigenvalues of the trap normal 
modes for different numbers of trapped ions A^. The results are shown in fig. (3). The 
function varies slowly with the value of A^, and , for N > 10, we can, to a good approx- 
imation, replace it by a constant S(A^) ^ 0.82. Thus we obtain the following upper 
limit on the total probability of the "extraneous" phonon modes becoming excited: 



-fp.r* ^ 



'2.6^0^' 



Nu 



(6.17) 



Thus we obtain the following sufficiency condition for the validity of Cirac and ZoUer's 
Hamiltonian 



'2.6^0^' 



Nu 



<1. 



(6.18) 



7 Conclusion 



In this preceeding sections we have reviewed the theoretical basis for cold trapped ion 
quantum computation. How these various laser-ion interaction effects may be combined 
to perform fundamental quantum logic gates is descibed in the seminal work of Cirac 
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and ZoUer [Q]. Using the formulae given here one can determine, for example, the laser 
field strength required or the separation between ions in the trap. Such things are of 
great importance in the engineering of practical devices. 

Finally there is the question of what type of ion to use. Figure 4 show the energy 
levels of four suitable species of ion. These have been chosen based two criteria: that 
the lowest excited state has a forbidden transition to the ground state, and their 
popularity amoungst published ion trapping experiments. It is not intended that this 
is an exhaustive list of suitable ions, but rather it is to show the properties of typical 
ions. 
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Appendix 

The vectors cf are usual normalized spherical basis vectors: 

c« = --^(l,-.,0), (A.l) 

c(°) = (0,0,1), (A.2) 

^ 1,^,0). (A.3) 



,(-1) 



V2 



(A.4) 



Note that 



c 



cW = (-l)9c(-^)* (A.5) 



The second rank tenors c\f are given by the formula 



y 



' V 3 ^ ^ £- _ \mi 1712 -q I ' ' 



15 



Explicity these five tensors are: 



Note that 





i 




1 




1 
1 

3 
\ 




1 




1 




«J 



(A.8) 

(A.9) 

(A.IO) 

(A.ll) 

(A.12) 



^ = (-1)^4'^^* (A.13) 

EcSlH^S'T = IV- (A.14) 
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Figure Captions 

Figure 1. A schematic diagram of ions in a linear trap, to illustrate the notation used 
in this paper. 



Figure 2. The relationship between the number of trapped ions A^ and the minimum 
separation. The curve is given by ( |2.8|) while the points come from the numerical 
solutions of the algebraic equations 



Figure 3. The function S(A^) defined by ( |6.16|) . 



Figure 4. Energy level diagrams for four species of ions suitable for quantum compu- 
tation. Wavelengths and lifetimes are given for the important transitions, the numbers 
in square brackets being the reference for the data. The lifetime is the reciprocal of the 
Einstein A coefficient defined in ( p.9| ) and ( p.lO| ). The thick lines are dipole allowed 



(El) transitions, the thin lines quadrupole allowed (E2) transitions. The atomic num- 
ber and the mass number of the most abundant isotope (with its relative abundance) 
are also given. Since all of these isotopes have an even mass number, they do not have 
a nuclear spin. 



Table Captions 

Table 1. Scaled equilibrium positions of the trapped ions, for different total numbers 
of ions. This data was obtained by numerical solutions of ( p.5|) . The length scale is 
given by 



Table 2. Numerically determined eigenvalues and eigenvectors of the matrix Anm de- 
fined by ( p.3|) , for 2 to 10 ions. The eigenvectors are normalized as defined by (|3.6|) . 
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i 



Mirror 



fY 




■^x 



Laser / 
Beam 



ions 



N Scaled equilibrium positions 



2 -0.62996 0.62996 

3 -1.0772 1.0772 

4 -1.4368 -0.45438 0.45438 1.4368 

5 -1.7429 -0.8221 0.8221 1.7429 

6 -2.0123 -1.1361 -0.36992 0.36992 1.1361 2.0123 

7 -2.2545 -1.4129 -0.68694 0.68694 1.4129 2.2545 

8 -2.4758 -1.6621 -0.96701 -0.31802 0.31802 0.96701 1.6621 2.4758 

9 -2.6803 -1.8897 -1.2195 -0.59958 0.59958 1.2195 1.8897 2.6803 

10 -2.8708 -2.10003 -1.4504 -0.85378 -0.2821 0.2821 0.85378 1.4504 2.10003 2.8708 



1.2 



03 
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10 



20 30 

Number of Ions, N 



40 



50 



N=2 



N=5 



Eigenvalue 


Eigenvector 






1 


(0.7071, 0.7071) 




3 


(-0.7071, 0.7071) 




N=3 








Eigenvalue 


Eigenvector 




1 


(0.5774, 0.5774, 0.5774) 




3 


(-0.7071, 0, 0.7071) 




5.8 


(0.4082, -0.8165, 0.4082) 




N=4 






Eigenvalue 


Eigenvector 


1 


(0.5, 0.5, 0.5, 0.5) 


3 


(-0.6742, -0.2132, 0.2132, 0.6742) 


5.81 


(0.5, -0.5, -0.5, 0.5) 


9.308 


(-0.2132, 0.6742, -0.6742, 0.2132) 



Eigenvalue 


Eigenvector 


1 


(0.4472, 0.4472, 0.4472, 0.4472, 0.4472) 


3 


(-0.6395, -0.3017, 0, 0.3017, 0.6395) 


5.818 


(0.5377, -0.2805, -0.5143, -0.2805, 0.5377) 


9.332 


(-0.3017, 0.6395, 0, -0.6395, 0.3017) 


13.47 


(0.1045, -0.4704, 0.7318, -0.4704, 0.1045) 



N=6 



Eigenvalue 


Eigenvector 


1 


(0.4082, 0.4082, 0.4082, 0.4082, 0.4082, 0.4082) 


3 


(-0.608, -0.3433, -0.1118, 0.1118, 0.3433, 0.608) 


5.824 


(-0.5531, 0.1332, 0.4199, 0.4199, 0.1332, -0.5531) 


9.352 


(0.3577, -0.5431, -0.2778, 0.2778, 0.5431, -0.3577) 


13.51 


(0.1655, -0.5618, 0.3963, 0.3963, -0.5618, 0.1655) 


18.27 


(-0.04902, 0.2954, -0.6406, 0.6406, -0.2954, 0.04902) 



N=7 



Eigenvalue 


Eigenvector 


1 


(0.378, 0.378, 0.378, 0.378, 0.378, 0.378, 0.378) 


3 


(-0.5801, -0.3636, -0.1768, 0, 0.1768, 0.3636, 0.5801) 


5.829 


(-0.5579, 0.031, 0.3213, 0.4111, 0.3213, 0.031, -0.5579) 


9.369 


(-0.3952, 0.445, 0.3818, 0, -0.3818, -0.445, 0.3952) 


13.55 


(-0.213, 0.5714, -0.1199, -0.4769, -0.1199, 0.5714, -0.213) 


18.32 


(0.08508, -0.4121, 0.5683, 0, -0.5683, 0.4121, -0.08508) 


23.66 


(0.02222, -0.1723, 0.4894, -0.6787, 0.4894, -0.1723, 0.02222) 



N=8 



Eigenvalue 


Eigenvector 




1 


(0.3536, 0.3536, 0.3536, 0.3536, 0.3536, 0.3536, 0.3536, 0.3536) 




3 


(-0.5556, -0.373, -0.217, -0.07137, 0.07137, 0.217, 0.373, 0.5556) 




5.834 


(-0.5571, -0.0425, 0.2362, 0.3634, 0.3634, 0.2362, -0.0425, -0.5571) 




9.383 


(0.4212, -0.3577, -0.4093, -0.1647, 0.1647, 0.4093, 0.3577, -0.4212) 




13.58 


(-0.2508, 0.5479, 0.0669, -0.364, -0.364, 0.0669, 0.5479, -0.2508) 




18.37 


(0.1176, -0.4732, 0.4123, 0.3039, -0.3039, -0.4123, 0.4732, -0.1176) 




23.73 


(-0.04169, 0.2703, -0.561, 0.3324, 0.3324, -0.561, 0.2703, -0.04169) 




29.63 


(-0.009806, 0.09504, -0.3398, 0.6127, -0.6127, 0.3398, -0.09504, 0.009806) 




N=9 




Eigenvalue 


Eigenvector 


1 


(0.3333, 0.3333, 0.3333, 0.3333, 0.3333, 0.3333, 0.3333, 0.3333, 0.3333) 


3 


(-0.5339, -0.3764, -0.2429, -0.1194, 0, 0.1194, 0.2429, 0.3764, 0.5339) 


5.838 


(-0.5532, -0.09692, 0.1658, 0.3078, 0.3531, 0.3078, 0.1658, -0.09692, -0.5532) 


9.396 


(-0.4394, 0.2828, 0.4019, 0.2558, 0, -0.2558, -0.4019, -0.2828, 0.4394) 


13.6 


(0.2812, -0.5108, -0.1873, 0.2228, 0.3881, 0.2228, -0.1873, -0.5108, 0.2812) 


18.41 


(0.1465, -0.5015, 0.2582, 0.4005, 0, -0.4005, -0.2582, 0.5015, -0.1465) 


23.79 


(0.06133, -0.3407, 0.5274, -0.02271, -0.4505, -0.02271, 0.5274, -0.3407, 0.06133) 


29.71 


(-0.01969, 0.1639, -0.4614, 0.5098, 0, -0.5098, 0.4614, -0.1639, 0.01969) 


36.16 


(-0.004234, 0.05021, -0.2195, 0.4939, -0.6408, 0.4939, -0.2195, 0.05021, -0.004234) 



N=10 



Eigenvalue 



Eigenvector 



1 



(0.3162, 0.3162, 0.3162, 0.3162, 0.3162, 0.3162, 0.3162, 0.3162, 0.3162, 0.3162) 



(-0.5146, -0.3764, -0.26, -0.153, -0.05056, 0.05056, 0.153, 0.26, 0.3764, 0.5146) 



5.841 



(-0.5476, -0.1382, 0.1079, 0.2544, 0.3235, 0.3235, 0.2544, 0.1079, -0.1382, -0.5476) 



9.408 



(0.4524, -0.2189, -0.3786, -0.3024, -0.1123, 0.1123, 0.3024, 0.3786, 0.2189, -0.4524) 



13.63 



(0.3059, -0.4689, -0.2629, 0.09726, 0.3287, 0.3287, 0.09726, -0.2629, -0.4689, 0.3059) 



18.45 



(0.1721, -0.5098, 0.1267, 0.3959, 0.194, -0.194, -0.3959, -0.1267, 0.5098, -0.1721) 



23.85 



(0.08046, -0.3902, 0.4528, 0.1795, -0.3225, -0.3225, 0.1795, 0.4528, -0.3902, 0.08046) 



29.79 



(0.03062, -0.2232, 0.505, -0.3078, -0.3154, 0.3154, 0.3078, -0.505, 0.2232, -0.03062) 



36.26 



(-0.009023, 0.09371, -0.338, 0.5419, -0.2886, -0.2886, 0.5419, -0.338, 0.09371, -0.009023) 



43.24 



(0.001795, -0.0256, 0.134, -0.3656, 0.5897, -0.5897, 0.3656, -0.134, 0.0256, -0.001795) 



^4 



0.8 



0.6 



0.4 



0.2 
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20 30 

Number of Ions, N 



40 



50 



Calcium II 

Atomic Number 20 
Mass number A = 40 (96.7%) 



4^P 



3/2 



4^Pi 



/2' 



\ 



854.209nm [18] 
101 nsec [21] 



396.847nm [18] 
7.7±0.2 nsec [20] 



849.802nm [18] 
901 nsec [21] 



866.214nm [18] 
94.3 nsec [21] 

\ 




3^D 



5/2 



3^D 



3/2 



QQQQ^^ riQi 729.147 nm 18 
393.366nm [18] . .. ..q. 

7.4±0.3 nsec [20] ^-^^sec. [i^J 




732.389 nm [18] 
16 sec. [19] 



Strontium II 

Atomic Number 38 
Mass number A = 88 (82.6%) 



5^P 



3/2 



5^Pi 



/2' 



\ 

1033.01nm [22] 
115 nsec [3] 



421.6706nm[22] 
7.87 nsec [24] 



1003.94nm [22] 
901 nsec [3] 



1091.79nm [22] 
105 nsec [24] 

\ 




407.886nm [22] 
6.99 nsec [24] 




4^D 



5/2 



7 

674.025589 nm [23] 
345±33 msec. [25" 



4^D 



3/2 



0066 nm [22] 
395±38 msec. [25] 



Barium II 

Atomic Number 56 
Mass Number 138 (71.7%) 



6^P 



3/2 



6^ Pi 



/2' 



493.4 nm [26] 
ll±lnsec[24] 



\ 

614.2nm [26] 
27±3 nsec [24] 



585.4 nm [26] 
210±30nsec[24] 

\ 



649.7 nm [26] 
30±4 nsec [24] 



\ 




455.4 nm [26] 
8.5±0.6 nsec [24] 



5'D 



5/2 



5^D 



3/2 



1.761 |im [26] 
47±16 sec. [27] 



2.051 |im [26] 
17.5±4.0 sec. [28] 



Mercury II 

Atomic Number 80 
Mass Number 202 (29.8%) 



{5di06p}'P3/2 



{Sd'HpfF,/ 



1/2 \ 



\ 



T 



991.4 nm [34] 
330±50 nsec [30] 



10.67 |im [34] 
3.0±0.5 iisec [30] 



{5d96s2}'D3/2 



398.0 nm [6] 
164.9 nm [29] 4.0±0.6 nsec [30] 
0.95±0.07 nsec [30] V 



|9/;.2i2 



{5d^6s^}^D5/2 



194.2 nm [29] 
2.3±0.3 nsec [30] 



197.8 nm [33] 
0.020±0.002sec. [31] 



{5di0 6p}'Si/2 




281.5766±0.00005nm[32] 
0.098±0.005sec. [31] 



